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Abstract 

We introduce the boundary conditions corresponding to the imaginary- 
time (Matsubara) formalism for the finite-temperature partition func- 
tion in d + 1 dimensions as constraints in the path integral for the 
vacuum amplitude (the zero-temperature partition function). We im- 
plement those constraints by using Lagrange multipliers, which are 
static fields, two of them associated to each physical degree of free- 
dom. After integrating out the original, physical fields, we obtain an 
effective representation for the partition function, depending only on 
the Lagrange multipliers. The resulting functional integral has the 
appealing property of involving only d-dimensional, time independent 
fields, looking like a non local version of the classical partition func- 
tion. We analyze the main properties of this novel representation for 
the partition function, developing the formalism within the context of 
two concrete examples: the real scalar and Dirac fields. 



1 Introduction 

About fifty years ago, in a pioneering work [T], the theoretical foundations 
for a systematic treatment of Quantum Field Theory (QFT) at finite tem- 
perature (T > 0), were laid down. That approach, the now called Matsubara 
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(or imaginary time) formalism has been very successful indeed in allowing 
for the evaluation of thermal effects in QFT, both in the High Energy 1 and 
Condensed Matter Physics realms. 

Another landmark in the subsequent development of that formalism was 
introduced in [3], where Matsubara's work was extended to relativistic QFT. 
The crucial discovery of the time-periodicity (antiperiodicity) conditions for 
the bosonic (fermionic) Greens's functions opened the door to a quite natural 
(and fruitful) extension of many concepts and ideas originally introduced at 
T = 0, to the finite temperature context. To name just a few examples: the 
thermal Ward-Takahashi relations, Goldstone's theorem, the Kubo-Martin- 
Schwinger (KMS) relations and the renormalization group at T > 0. Besides, 
the notion of Abelian and non Abelian gauge fields, with all its consequences 
for particle physics [3 [HI could also be set up in the finite temperature 
context in a quite natural way. 

In the path-integral framework for the imaginary time formalism, ther- 
mal effects depend strongly on the fact that the imaginary-time coordinate 
is compactified. That is, fields corresponding to a d + 1-dimensional theory, 
must be defined on S 1 x R d , where the radius of S 1 is proportional to the 
inverse temperature. This fact has, of course, many important consequences 
for the construction of topologically non-trivial solutions jH| since they nat- 
urally depend upon the global properties of the spacetime manifold (unlike, 
for example, the UV structure of the theory). On the other hand, the topo- 
logically trivial sector is also affected by the fact that the frequencies become 
discrete variables (the Matsubara frequencies). 

A feature of the Matsubara formalism (shared by the real-time formula- 
tion) is that a time dependence for the fields is unavoidable, even when one 
limits oneself to the calculation of static thermal properties. With the aim 
of constructing a new representation where only static fields are involved, in 
this paper, we introduce an alternative way of dealing with finite-temperature 
QFT. It is inspired by a recent paper in which a constrained functional in- 
tegral approach is used to implement the effect of fluctuating boundaries 
in the Casimir effect 0. In the present context, this allows one to intro- 
duce the periodicity conditions by means of Lagrange multipliers (which are 
d-dimensional when the field lives in d + 1 dimensions). Then the origi- 
nal fields can be integrated, what leaves a functional depending only on the 
(i-dimensional Lagrange multipliers. Besides the economy due to the 'dimen- 
sional reduction', there are also some new interesting properties, as we shall 
show in the different sections of this article. 

The structure of this paper is as follows: in section we introduce the 

^ee, for example 
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constraint formalism to deal with T > 0, for the simplest possible example, 
i.e., the harmonic oscillator. Equipped with these results, the case of the 
real scalar field is dealt with in section 03 both for the free and the self- 
interacting cases. In section HI we consider the Dirac field, and in El we 
present our conclusions. 



2 Harmonic oscillator 

We shall introduce the method by considering the simplest (yet non triv- 
ial) example, the harmonic oscillator. In the subsequent sections, the same 
approach shall be extended, in a quite natural way, to interacting QFT's. 

Let us begin by writing the phase-space path integral [E] for Z , the 
vacuum functional for the harmonic oscillator: 

Z = J VpVqexp{-So[q(T),p(r)]} , (1) 
where So is the imaginary time first-order action: 

S [q(T),p(r)} = / dT[-ip(T)q(T)+Ho(p(T),q(T))]\ (2) 



with Ho denoting the classical Hamiltonian, 

n (p,q) = \{p 2 +u;V) (3) 

(for simplicity, we have assumed that the oscillator has unit mass). We have 
introduced in (|T} the formal expression for the phase-space path integral 
integration measure: 

VpVqs n (4 ) 

— 00<T<00 

whose precise definition may, of course, be given in terms of a discrete ap- 
proximation to the path integral 2 . 

On the other hand, in the Matsubara formalism, Z (/3), the thermal par- 
tition function (for the same system) can be written as: 

Zo{(3) = ( VpVqe^ dT ^-W^ 2 ^\ , (5) 

J periodic 



2 The reason why the phase space path integral (rather than the more usual configura- 
tion space one) is used, will become clear at the end of this section. 



3 



where the integration is taken over periodic paths q(f3) = q(0), p(/3) = p(0) 3 
and we have used units such that (3 — ^, i.e., ks = 1. 

To simplify matters, we shall measure the energies with respect to the vac- 
uum (zero energy). For the partition function, this means that Zq{oo) = 1; 
hence, we shall absorb temperature-independent factors by introducing a 
suitably chosen normalization constant in Z (f3). 

To reproduce 2 (/3), we shall follow a quite straightforward procedure: 
we will impose the periodicity conditions for the phase space variables as 
constraints for the integration variables in Z . To that end, we write: 

Zo(j3) = N I VpVq 5(q{{3) - q(0)) 5{p{(3) -p(0)) e - So[q{T) ^ r)] , (6) 



where Af is a constant, introduced to enforce the condition Z (oo) = 1, and 
Sq is the action defined in (|2j). 

The imaginary-time r runs from — oo to +oo, while the constraints iden- 
tify just the fields at and (3. Thus one should expect the emergence of 
a ^(oo) factor in (JUJ), due to the integration of fields outside of the inter- 
val [0, 0\. That factor is, however, /3-independent, and shall therefore be 
absorbed by N . 

To proceed, we exponentiate the two 5-functions, by means of (two) La- 
grange multipliers £1 and £ 2 (they are just real variables in this case); we are 
thus lead to a new expression for Zo((3), which may be written as follows: 

if + ^dr ,(t)€i (<5(t-/3)-5(t)) +p(r) & (<5(r-/3)-5(r)) 

x e L '\ 



The integral over p and q (which we shall perform first) is obviously Gaus- 
sian. Indeed, interchanging the order of integration, we see that JZj) may be 
rewritten as follows: 

where we have introduced Q = (Q a ) (a = 1, 2) for the canonical coordinates, 

such that Qi = q and Q 2 = p. Also, j a ( T ) = £a (£( r — 0) — ^( r )), an d KL^ 

3 This symmetric form of the boundary conditions in the path integral for Zq{(5) is 
explained in 
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are the elements of the 2x2 operator matrix /C, given by: 

The result of integrating out Q a , may be written as follows: 
Z Q ((3) = AT27r(det£)^ J 0^ e ~¥^ h ^ 



(10) 



with 

M = 2fi(0) - fi(/3) - n(-/3) , (11) 

where fi(r) denotes the inverse of the operator /C of (jOJ); its explicit form 
may be easily found to be: 



v ; V -fsgn(r) f J K J 

(sgn = sign function). Equation (fT2|) can be used in ([TTf . to write: 

M = ( W ~ 1 I ) (n B (o;) + I)" 1 , (13) 

where 

n B (cu) = (e^ - l)- 1 (14) 

is the Bose-Einstein distribution function. 

The resulting integral over the £ a variables becomes: 

Z B = —e ^B^+n , (15) 

J 27T 

where we have set: 

J\f = (det£) 5 , (16) 

by the normalization convention. The resulting integral is over the two (real) 
variables £ a , which are 0-dimensional fields, one dimension less than the + 1 
dimensional original theory. 

It is interesting to see how this integral may be compared with the one 
corresponding to classical statistical mechanics. To that end, we evaluate the 
partition function in the classical (high-temperature) limit. In that limit we 
approximate the integrand and Z (f3) becomes: 

mp) - / Cr e "^ (Sl ' 6) o*«i). (17) 



where: 

#(6,6) = + (is) 

We see that (fTTj) corresponds exactly to the classical partition function for a 
harmonic oscillator, when the identifications: £1 = p (classical momentum), 
and £2 = <? (classical coordinate) are made 



UP) 




On the other hand, had the exact form of the integral been kept (no ap- 
proximation), we could still have written an expression similar to the classical 
partition function, albeit with an 'effective Hamiltonian' -ffe//(£i, £2): 

UP) = I ^e-^ff^, (20) 

where: 

This shows that the quantum partition function may also be written as a 
classical one, by using a /3-dependent Hamiltonian, which of course tends to 
its classical counterpart in the high-temperature limit. This representation 
is also valid for interacting theories, as shown in 13.21 

By integrating out the Lagrange multipliers in the (exact) expression for 
the partition function fTHj) . we obtain: 

Z Q (P) = n B (u) + 1 = - _ ^ . (22) 

which is the correct result. 

We emphasize an important fact that has emerged from an analysis of the 
classical (high-temperature) limit in this case, namely, the Lagrange multi- 
pliers have a physical interpretation in the classical limit: the one associated 
to the periodicity condition for q plays the role of the classical momentum, 
while the one corresponding to the periodicity for the momentum becomes 
a generalization of the classical coordinate. The same interpretation might 
also be retained far from the classical limit, but then the Hamiltonian departs 
from the classical one, receiving quantum corrections. 

To conclude this section, we mention the fact that the (apparently) sim- 
pler procedure of working in terms of the configuration space path-integral 
yields an incorrect answer. Indeed, had one used that formulation, one should 
have needed to impose the periodicity constraint just for q(r): 

Z {P) = Af J Vq5(q(f3)-q(0)) x e - * ^ <^0 2 +"V) . (23) 
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The 5-function could then have been easily exponentiated by using just one 
Lagrange multiplier, X'- 



Zo(/3) = N I Vq I ^e^Me-^-^M . 
Performing the Gaussian integral over q(r) would have lead to: 

/+oo 
d X e~ Q( *> , (25) 
-oo 

where: 

Q(X) = u~ l [n B (u) + l]~ l x 2 . (26) 
Choosing Af~ l = Z (oo), and integrating out x> we would have seen that 

MP) = „ 1 g (wrong) (27) 

which is of course incorrect, due to the power 1/2 (we have neglected (3- 
independent factors). The reason why this formulation yields an incorrect 
answer may be traced back to the fact that it corresponds to a different con- 
straint. Indeed, the previous, phase space treatment, amounts to enforcing 
periodicity for the coordinates and their momenta (essentially, their time 
derivatives). Thus, it corresponds to requiring the periodicity for q(r), and 
also the time independence of that constraint. 



3 Real scalar field 



3.1 Free partition function 

In what is perhaps the most straightforward extension to QFT of the previous 
section's results, we shall consider here a real scalar field (p in d+1 (Euclidean) 
dimensions: f(x) = (f(r, x) where x = (r, x) G M^ +1 ^, r G K and x G R". 

The free Euclidean action, So, defined in terms of the phase-space vari- 
ables is 

' d d+l x 



So 



- md T tp + H (ir, ip) 



with 



H (7T,^) EE ^7T 2 + |V<^| 2 + mV 



(28) 



(29) 
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We then have to implement the periodic boundary conditions both for 
<p(T, x) and its canonical momentum vr(r, x) 



<p(P,x) = V?(0,x) , Vx G 
tt(/3,x) = tt(0,x) , Vx G , (30) 

which requires the introduction of two Lagrange multiplier fields, as in the 
previous section about the harmonic oscillator. We denote them by £ a ( x ), 
a = 1, 2 (they are r-independent). 

To simplify matters, we now define a two-component field <3> = ($ a ), 
a = 1, 2, such that $i = (p and $2 = 7r. Then, an entirely analogous proce- 
dure to the one introduced for the harmonic oscillator yields: 



(31) 



where: 

j a (x) = £ a (x) (S(t - (3) - S(t)) . (32) 
Note that the measure for £ is formally given by: 

^nffwx^ (33) 

while: 

2tt 

r,x 

and the operator matrix K, is: 



tt d(p(r, x)d7r(r, x) . , s 

£>$ = J ^ v - ' — ^-J. , (34) 
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where we have introduced /i = V— V 2 + m 2 , the first-quantized energy oper- 
ator for massive scalar particles. 

Performing the integral over $, and fixing Af to verify the normalization 
condition, yields the partition function in terms of the Lagrange multipliers: 

with: 

M = 2fi(0) - £l(P) - tt(-p) , (37) 



and 



Then, 



where 



^ = ( v n {fiB + ' (39) 

n B = — • (40) 

e^-1 V 7 

Coming back to the expression for Zq{0), we see that: 

Z (/3) = Jv^exp[- 1 -Jd d xJ d d y[U*) (x^ 1 (n B + l)" 1 ^) £x(y) 
+ 6(x) (x|/ i (f iB + l)- 1 |y)6(y)]} • (41) 
By a simple field redefinition, we see that: 

Z (/3) = det (n B + 1) (42) 
which can be evaluated in the basis of eigenstates of momentum to yield: 

Zo(P) = II [M**) + 1] (43) 



where = -y/k 2 + m 2 . The free-energy density, Fq(/3), is of course: 

For the sake of completeness, we derive here the explicit form of the path 
integral for the classical, high-temperature limit. When (3 ~ 0, the integral 
for the partition function becomes: 

Z (P) ~ Jvte-? 11 ® , (45) 

where: 

H(0 = \j ^K(x) + |V6(x)| 2 + m 2 £ 2 2 (x)] . (46) 

This is, again, the usual classical expression for the partition function, with 
the Lagrange multipliers playing the role of phase space variables, and the in- 
tegration measure being the corresponding Liouville measure. Besides, there 
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is another fact that becomes here more transparent than in the harmonic 
oscillator case, and we would like to point out. It is that the representation 
ijHj) always involves static fields, unlike in the Matsubara formalism. The 
price to pay for this 'dimensional reduction' is that the resulting 'action' (the 
exponent of the functional to be integrated) is spatially non local. It becomes 
local only in the high-temperature limit. 



3.2 Self-interacting real scalar field 

When field self-interactions are included, instead of the free action So, we 
must consider instead 

S = S + Sj, (47) 

where the free action So, has already been defined in while Si is a 

self-interaction term. We shall assume it to be of the type: 

Si = / d d+1 xV(<p) , (48) 



V(ip) being an even polynomial in ip, with only one (trivial) minimum. Pro- 
ceeding along similar lines to the ones followed for the free field case in the 
preceding section, the partition function for the interacting system can be 
written in the form: 



(49) 



where $, as well as the 'current' j a have already been defined for the free 
case, in the previous subsection. The constant Af is again introduced to 
satisfy Z{po) = 1. On the other hand, since the fields are assumed to tend 
to zero at infinity, (3 — > oo implies that the term involving j vanishes in this 
limit. This means that 



Af' 



J Vi Jv<$> e- 3 ^ . (50) 



There are many different paths one could follow from now on in order to 
evaluate the partition function. We choose to adopt a procedure that makes 
contact with quantities defined for QFT at T = 0, in such a way that the 
T^O theory is built 'on top of it'. 

Indeed, recalling the definition of the generating functional for connected 
correlation functions, W, we may write: 



Af J V§ exp -S($) + % J d d+1 xj a <f> c 



= e~ w ^, (51) 
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so that the partition function Z((3) becomes: 



Z(j3f) = J Vi e~ w ^ . (52) 

We use the small j to denote the 2-component current which is a function of 
the Lagrange multipliers, as defined in (|32j). A capital J shall be reserved to 
denote a completely arbitrary 2-component source, so that: 



J\f J £>$ exp -£($) + i J d d+l xJ a $ c 



-W(J) 



(53) 



Note that this generating functional yields correlation functions includ- 
ing (f or 7r legs, so it would seem to contain much more information than 
the standard one. That is not so, however, since W( Ji, J-i) is related to 
W(Ji) = W(Ji,0), the standard generating functional of T = y?-field cor- 
relation functions: 

W(Ji) = -ln{ J V p e -fd d M^) 2 H™ 2 * 2 +v(v)-iJiv]y (54) 

That relation that may be obtained, for example, by integrating out the 
canonical momentum field tt (a Gaussian integral): 

1 ' ld+1, 



W(Ji, J 2 ) = W(Ji - z9 r J 2 ) + - J d d+l xJ 2 2 (x) . (55) 

The outcome of the previous relation is that correlation functions including 
7r legs may be obtained by combining correlation functions where all the 
canonical momenta have been replaced by scalar fields (at the same spacetime 
points), and derivatives of them. 

Besides, taking into account the fact that our normalization for Z{j3) 
implies that W(0, 0) = 0, we are ready to write a formula for the path- 
integral representation we were looking for. Indeed, defining H e ff(£), the 
'effective Hamiltonian' for £, by means of the expression: 

HefAO = iwtf), (56) 

we see that the partition function is given by: 

Z((3) = JVC exp \-(3H eff {0] ■ (57) 

This is one of the most important results in this article, since it yields the path 
integral for the quantum partition function as a classical-looking functional 



11 



integral, involving an effective Hamiltonian which takes into account all the 
T = quantum effects. Indeed, the usual functional expansion for W(J) is: 



W(J) 



oo 1 

^ n! 

n=2 



■ • ■ d 



d+1 



•^7 



{51 



where is i n times the n-point connected correlation function for the 

field $. Thus, knowing those correlation functions, at a given order in a loop 
expansion, one may obtain an analogous functional expansion for the effective 
Hamiltonian. Indeed, to do that, one should perform the integrations over 
the imaginary-time variables (taking advantage of the 5- functions): 



oo 

HefM) = £ ~ / ^ ■ ■ ■ dd *n HM^fa, . . . , X n ) 
n=2 ' 



x ^ 1 (xi)---^ B (x n ) 



(59) 



where 



H (n) 



\ J dr x ...dr n W£ >an (xi, x n ) 



x (5(n - P) - 6(n)) ... (8(r n - 13) - 5(r n )) . 

Of course, one will usually evaluate just a few of the terms in the ex- 
pansion for W, and as a consequence will obtain an approximation to the 
effective Hamiltonian. As an exercise, we evaluate the partition function that 
results from the first non-trivial approximation, namely, to keep only the first 
term, the 2-point correlation function 4 : 



where 



and 



^ J dxidx 2 ^o(3Ci)fl'S ) (xi l x 2 )£ 6 (x 2 ) 



(60) 
(61) 



1 r 



wi?(Ax i; 0,x 2 )-Wi?(0,x i; Ax 2 



(62) 



4 The meaning of this approximation shall be discussed at the end of this calculation. 
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Introducing Fourier transforms for the correlation functions W a & and the 
kernel, , we see that (f62l) implies: 

^(w,k) = -g y ^ — [1 - oos^wg^w.k). (63) 

(2) 

On the other hand, the kernels W ah can ^ e obtained from the T = 
correlation function W^y which is, essentially, the full quantum propagator 
for the y?-field: 

W[l\x 1 ,x 2 ) = - (ip{x x )<p{x 2 )) . (64) 

(2) 

Indeed, the explicit relations between the correlation functions W ab with 
a 7^ 1 or b ^ 1 and Wy , in momentum space, are: 

(u,,k) = 1 - j^WyVk) 
W{?^,k) = -ujW[f(cv,k) 

Wg } (o;,k) = coWg^.k)). (65) 

At this point, in order to evaluate the frequency integrals, we need of course 
to make some assumptions about the structure of w[f . We shall assume 
that there is only one stable particle, and that the renormalization condi- 
tions corresponding to the physical mass and wavefunction renormalization 

' (2) 

have been imposed. Then Wn has only two single poles, at the points 
±iE ren (k), where E ren (k) is the dispersion relation, with a renormalized 
mass, and including all the quantum corrections. 

Considering each component of equation (JHSj) for specific values for a and 
6, we see that their respective to integrals are easily evaluated under the 
previous assumptions. The results are: 

(66) 



H 



(2), 
11 I 



E^ n (k) n B (E ren (k)) 



Hf 2 \k) = E ren (k) n B (E ren (k)) 



(67) 



while H 



(2) 
12 



H. 



(2) 
21 



0. 



Since the functional integral has the same structure as the one for the 
free case, to evaluate Z(f3) in this approximation we just need to modify the 
expression for the energy in the free field result: 



(68) 
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where 

Z^\(3) = J] [n B (E ren (k)) + 1] (69) 

k 

which yields the free energy 



(70) 



Of course, this expression looks like a the free energy for a free field, but 
with the renormalized energies instead of the free ones. This approximation 
amounts, of course, to considering the particle states as non interacting, after 
the self-energy corrections have been taken into account (in the 2-point, T = 
correlation function). A possible way to understand this approximation is 
by the large- N limit, whereby the correlation functions involving more than 
2 points are suppressed, and the approximation used here becomes exact. 

An interesting approximation scheme, related to the functional expansion, 
is given by the High-Temperature expansion. In terms of H e jj, we see that 
the only place where (3 appears is in the source j. Performing an expansion 
of j in powers of j3: 

00 on 

Ja(x) = £ a (x) E^ 1 )"^"^)' ( 71 ) 
n=0 

and inserting this expansion into the expression for H e ff, we see that the 
kernels corresponding to its functional expansion become: 

°° /?niH hnm-1 

fx, ... x 1 = V - 

ni—l,— ,n m 

x [d T1 ---d Tm W^. a Jx 1 ,--- ,x m )] T ^ . (72) 

We see from this expression that, when (3 — > 0, the leading term is linear 
in j3 and only survives: 

^) 2 ( Xl ,X 2 ) ~ P 2 [d T AM%Mi^)] T1 ,r^0 ■ ( 73 ) 

We conclude this section mentioning the use of higher order approxima- 
tions to the effective Hamiltonian, say, by keeping N terms in the expansion: 

N 

n=2 ' J 

x , (74) 
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(where only even values on n are non- vanishing in the sum). The result- 
ing partition function naturally admits a perturbative expansion, with the 

(2) 

quadratic kernel H ab determining the propagator, and the higher order terms 
defining nonlocal vertices. 



4 Dirac field 

We shall derive here the effective Hamiltonian for the partition function cor- 
responding to a massive Dirac field in d + 1 spacetime dimensions. The 
procedure will be essentially the same as for the real scalar field, once the 
relevant kinematical differences are taken into account. 

As in the previous section, we shall first deal with the free field case. 



4.1 Free Theory 

The action Sq for a free Dirac field ip in d + 1 Euclidean dimensions is: 

= J d d+1 x4){x){ft + m)4){x) , (75) 

where ft = 7 M <9 M , 7+ = 7 M and {7^,7^} = 25^. 

The corresponding T = partition function, Zo (/?)|/3->oo, is then given 

by: 

z£(oo) = J Vrj)V$e- s °®M, (76) 

where ip,ip are Grassmann-valued spinorial fields. 

To derive the thermal partition function, Zq(/3), we must impose antiperi- 
odicity conditions for both fields: 

tf>(/3,x) = -V(0,x) 

V(/5,x) = -</>(0,x) (77) 

as constraints on the Grassmann fields in the path integral (f76|) . Those 
conditions lead to the introduction of the two 5— functions: 

x exp[-3f(^)] . (78) 

Note that, since the Dirac action is of the first-order, the introduction of two 
constraints, and two Lagrange multipliers, appears in an even more natural 
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way than for the previous case. Those auxiliary fields, denoted by x( x ) and 
x(x) must be Grassmann spinors depending on the spatial coordinates only. 
The resulting expression for Zq(@) can then be written as follows: 



Z£(/3) = M / WWe 



(79) 



where r\ and 77 are (Grassmann) sources depending on x an d X through the 
relations: 



rj(x) = x (x) [S(t - (3) + *(t)] 
jfts) = x(x) [6(t - (3) + *(r)] 



(80) 



and the constant M is determined by imposing the usual normalization con- 
dition: Zq (oo) = 1. 
Integrating out ip,ip: 



VxVx exp 



/ d a+1 yfj(x){x\{0 + m)- l \y)rj(y) 



ld+1 



where we have set Af 



VxVx exp 
Zq (oo) and 



P H eff(x,X) 



(81) 



H eff (x,x) = j d d x J Ax(x)^ (2) (x,y)x(y) (82) 



with: 



tf (2) (x,y) = (x,0|(^ + m)^ 1 |y,0) + (x,/3|(^ + m)- 1 |y,/5) 
+ (x,0|(^ + m)- 1 |y,/3) + (x,/3|(^ + m)- 1 |y,0) 

= ±[2S f (0, x - y) + S f (/3, x - y) + S f ( - /3, x - y)] (83) 

In the last line, S/, denotes the Dirac propagator. Then a quite straightfor- 
ward calculation shows that 



#(x,y) = -(x\u(l-n F ) x |y > 



where 



1 + e 



j3!i 



(84) 



(85) 
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is the Fermi-Dirac distribution function, written in terms of h, the energy 
operator (defined identically to its real scalar field counterpart), u is a unitary 
operator, defined as 



. _ h D 
u = , h D = 7 ■ V + m . 
h 



Then we easily verify that: 

Z*(P) = detu det _1 [(l -h F )l] , 
(I = identity matrix in the representation of Dirac's algebra) 



4W 



-f3E(p) 



(86) 



(87) 



(88) 



with E(p) = -\/p 2 + m 2 and r d = dimension of the representation. 
We have used the fact that the determinant of u equals unity: 



det u 



J det 



l l ' P + 111 



p2 _|_ m 2j 



(89) 



Again, the procedure has produced the right result for the partition func- 



tion. 



4.2 Dirac field in a static external background 



We begin by highlighting the main differences which appear when an external, 
static (we deal with equilibrium thermal QFT), Abelian gauge field is coupled 
to the Dirac field. The first is that the action becomes: 



S f A) 



d 



d+l 



X 



ifj(x)( 



icy ■ A(x) 



m 



X 



(90) 



where we have used the A = gauge. The assumed r— independence of A, 
together with our gauge choice, allows us to carry on the derivation described 
for the free case, with minor changes, arriving to the expression: 



Z f ((3) = detu(A)det- 1 (n F (A)l) 
= e^det [(l + e-^ A ')l 



(91) 



where 



h(A) = 



77T 



D ee V - ieA 



(92) 
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and: 



det (7 • D + m) 
det V-D 2 + m 2 



(93) 



Notice that the factor det 



can be formally diagonalized in 



(l + e-W)l 

terms of the energies E\(A) in the presence of the external field. Thus we 
arrive to the expression: 



e iK(A) x 



-I3E X (A) 



(94) 



The factor e lK ( A > , on the other hand, is topological in origin, as it depends 
on the phase, K(A), of the determinant of hp, an operator that may be 
regarded as a Dirac operator in one fewer dimension. For Dirac fermions, 
we know that the phase of det hp can be non-trivial only when d is odd, 
i.e., when d + 1 is even. However, the 7-matrices appearing in det hp form 
a reducible representation of the Dirac algebra in d dimensions, with the 
matrix 7 T relating every eigenvalue to its complex conjugate. Thus, as a 
result, the phase K(A) vanishes. Of course, a non-vanishing result may be 
obtained for other fermionic systems, like Weyl fermions for d + 1 = even. 
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5 Conclusions 



We have shown that, by introducing the 'thermal' boundary conditions as 
constraints in the Euclidean path integral for the vacuum functional, we may 
obtain a novel representation for the partition function. This representation 
may be thought of as an integral over the phase space variables, weighted by a 
Boltzmann factor corresponding to an effective quantum Hamiltonian, H e ff, 
which reduces to the classical one in the corresponding (high-temperature) 
limit. 

We analyzed the main properties of this representation for the cases of 
the real scalar and Dirac fields, two typical examples that have been chosen 
for the sake of simplicity. It is not difficult to generalize the representation 
to the case of systems containing fermions interacting with bosons. For 
example, assuming that Sty, ip, $) is the first order action corresponding to 
a real scalar interacting with a Dirac field, we define the T = generating 
functional W(C, (, J) by: 

W(C, (,J) = M J V^ThpThf) e~ si ^ )+i * dd+lx (c^C+J^a) ? ( 95 ) 

where the sources are arbitrary. Then, H e ff can be obtained from the ex- 
pression: 

Heff(ZX,Z) = jjW(v,V,j) , (96) 

where r), fj and J are (the already defined) functions of the Lagrange multi- 
plier fields x, X an d £■ 

We conclude by mentioning that r-dependent correlation functions can 
also be calculated, with minor changes in the formalism. Indeed, they may 
be obtained by adding an arbitrary time-dependent source to the current 
that depends on the constraints [TU] . 
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